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Abstract 

Phase separation has been observed by several different experiments and it is beheved to be 
closely related with the physics of cuprates but its exactly role is not yet well known. We propose 
that the onset of pseudogap phenomenon or the upper pseudogap temperature T* has its origin in 
a spontaneous phase separation transition at the temperature Tpg = T* . In order to perform quan- 
titative calculations, we use a Cahn-Hilliard (CH) differential equation originally proposed to the 
studies of alloys and on a spinodal decomposition mechanism. Solving numerically the CH equa- 
tion it is possible to follow the time evolution of a coarse-grained order parameter which satisfies 
a Ginzburg-Landau free-energy functional commonly used to model superconductors. In this ap- 
proach, we follow the process of charge segregation into two main equilibrium hole density branches 
and the energy gap normally attributed to the upper pseudogap arises as the free-energy potential 
barrier between these two equilibrium densities below Tpg. This simulation provides quantitative 
results in agreement with the observed stripe and granular pattern of segregation. Furthermore, 
with a Bogoliubov-deGennes (BdG) local superconducting critical temperature calculation for the 
lower pseudogap or the onset of local superconductivity, it yields novel interpretation of several 
non-conventional measurements on cuprates. 

PACS numbers: 74.72.-h, 74.80.-g, 74.20.De, 02.70.Bf 
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I. INTRODUCTION 
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The existence of the pseudogap in all family of high-temperature superconductors (HTSC 
has been verified by several different experimental techniques as discussed by many reviews 
Oj. As a consequence of many years of scientific effort, there is a solid consensus of its 
existence at least in the underdoped regime. On the other hand, there is currently no 
agreement on such basic facts as to its nature and origin. After its discoveryR Q, it 
was realized that some experiments detected the pseudogap temperature T* at very high 
values while others would place it just above the critical temperature T^. This is probably 
because different probes are able to detect different properties but, the fact is that this 
large discrepancy triggered a variety of different proposals. Just to mention a few ideas and 
works; Emery et al. jj] called the high T* as T^*, the crossover temperature at which charge 
inhomogeneities become well defined and the low T* as T2* and associated it with a spin 
gap and they both merged into at the slightly overdoped region of the phase diagram. In 
their review Timusk and StattQ] also presented a similar phase diagram but they related the 
lower pseudogap temperature to T* and the upper one also to a crossover temperature T^. 
The lower and the higher T* were also considered as the opening of a spin and a charge ga^ 
respectively [5|l. The lower T* was also attributed to superconducting phase fiuctuationsLSl 

nnnnri 

and many different experiments claimed to have detected such fiuctuations)7l. 181 191. IIOI. 
Thus the existence of the two pseudogaps in the cuprates has been compiled by several 
works as the result of many different data. In fact, analyzing the data from angle- 

resolved photoemission (ARPES) and angle-integrated photoemission (AIPES), Ino et al. 13| 
could distinguish not two but three different energy scales. 

Another controversial point is whether the pseudogap and the superconducting gap have 
the same origin or not. Tunneling spectroscopy Q, Q, ^| seems to show that the gap 
evolves continuously from the superconducting into the normal phase without any anomaly, 
suggesting that the pseudogap and superconducting gaps have the same origin. The com- 
mon origin was also supported by some ARPESll7| and scanning tunneling spectroscopy 
(STM)|l8| data. Muon spin rotating experiments |l9|] characterized T* as the pair formation 
line in agreement with the fluctuation theories of pre-formed superconducting pairs 0,0, 2ol |. 
These STM and ARPES experiments have also measured the pseudogap in the overdoped 
region in opposition to many others [l|, Q, 13] which the pseudogap temperature line appears 
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to fall a little beyond the optimum doping value. On the other hand, intrinsic (c-axis inter- 
plane) tunneling spectroscopy 2l|, UA, UM results against a superconducting origin of 
the pseudogap what was also confirmed by the same type of experiment in high magnetic 
field |22|. This conclusion, against the common origin of the pseudogap and superconduct- 
ing gap, is also shared by Tallon and Loram after the analysis of data from many different 
experiments yj]. 

The above resumed paragraphs intended to show that, despite the enormous experimental 
effort after all these years, there are still some basic open questions in this field. These open 
questions motivated us to make the present novel work which connects the large pseudogap 
T* to the onset of phase separation. There is now considerable evidence that the tendency 
toward phase separation or intrinsic hole clustering formation is an universal feature of doped 



cuprates 
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271 128| . Phase separation in hole rich and hole poor regions was theo- 
and has been observed in the form of stripes j30|, ISlI and in the form 



retically predicted 
of microscopic grains or mesoscopic segregation by STM measurements 



Although 



the STM results has been questioned as a surface phenomena which does not reflect the 



nature of the 
diffraction 
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k electronic state[34], the inhomogeneities has also been seen by neutron 
, [3^ which is essentially a bulk type probe in underdoped and optimally 
doped region of the La2-xSrxCuOi phase diagram. Another bulk-type measurement using 
nuclear quadrupole resonance (NQR) 36| has observed an increase in the hole density spatial 
variation of La2-xSrxCu04 compounds (with 0.04 < a; < 0.15) as function of the tempera- 
ture. Despite these evidences, the majority of the theoretical approaches are based on the 
assumption that the holes are homogeneously doped into CuO planes, probably due to the 
argument that, in principle, macroscopic phase separation is prevented by the large Coulomb 
energy cost of concentrating doped holes into small regions. On the other hand, the above 
cited references are just a few of the large number of works which have detected some type 
of inhomogeneities in cuprates which seems to be intrinsic since it is present even in the best 
single crystals j^]. There are also experimental evidences for an intrinsic phase separation 
and cluster formation in many other materials like, for instance, manganites which are be- 



lieved to be another strong correlated electron materials 



38L |39[ and on rutheno-cuprates 



superconductors j40|. In fact, it has been argued that phase separation might be stronger in 



manganites 



33] than in cuprates. 



In this article we develop a novel approach to this issue as we apply to the large pseu- 
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dogap T* the theory of phase-ordering dynamics, that is, the growth of domain coarsenin 



when a system is quenched from the homogeneous phase into an inhomogeneous phase |41|. 
This phenomenon is also known as spinodal decomposition. One of the leading models de- 
vised for the theoretical study of this phenomenon for a conservative order parameter is 



based on the Cahn-Hilliard formulation 



42j |. The Cahn-Hilliard (CH) theory was originally 



proposed to model the quenching of binary alloys through the critical temperature but it 
has subsequently been adop ted to model many other physical systems which go through a 
similar phase separation|4l[ 13. 143^ . We show how the CH equation is derived from a typical 
Ginzburg-Landau (GL) free energy for a typical (conserved) order parameter, which is easily 
related with the density of holes, using an equation for the conservation of the order param- 
eter current. The CH equation is solved numerically by adopting a very efficient method 
(compared with usual first order Euler methods) semi-implicit (in time) finite difference 



scheme proposed by Eyre 



43l |. The numerical details have been analyzed elsewhere 



The main purpose to solve the CH equation for the hole density field and take the large 
pseudogap temperature T* as the phase separation temperature Tpg is that we can make 
quantitative calculations and get some insights on various HTSC non-conventional features: 
as the temperature goes down below Tp^, the distribution of hole density for a given com- 
pound evolves smoothly from an initially random variation taken as a Gaussian distribution 
around an average density p, since a purely uniform distribution does not segregate, into a 
kind of bimodal distribution. These simulations are used to demonstrate the charge inhomo- 
geneity and the stripe pattern formation in a square lattice as shown below. The pseudogap 
energy Eg or the large pseudogap temperature T* arises naturally as the GL potential barrier 
between the two equilibrium density phases, changes smoothly as the temperature decreases 
and reaches the maximum phase separation near zero temperature. If Tpg vanishes at a 
critical average means that all 

the compounds with average p < Pc may undergo a phase separation and evolves contin- 
uously into a complete separation characterized by a bimodal distribution with two major 
equilibrium densities (p+ and p-). For underdoped samples the phase separation is more 
pronounced, since Tpg is very large for these compounds. The difference between p^ and p_ 
should decrease for compounds with increasing average hole density p and the sharp peaks 
evolves into rounded peaks near p^. This provides an explanation for the neutron diffrac- 
tion data on the Cu — O bond length distribution and the observation of charge and 
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spin separation into stripe phases. On the other hand, the increase of the inhomogeneity 
(variation in p) as the temperature is decreased for a given sample was observed by the 
NQR experiments 36|], in agreement with the CH theory of the spinodal decomposition. On 
the other hand, these local differences in the charge distribution generate local microscopic 
(or mesoscopic) regions with different superconducting transition temperatures. The onset 
of local superconductivity may be identified as the lower pseudogap temperature or the 
temperature where the superconducting pairs start to appears. This second pseudogap has 



also been interpreted as the mean field temperature T by Emery and Kivelsonj iM. As 
the temperatures goes down between this lower T* and Tc more superconducting regions or 
superconducting droplets appear, they grow in size and quantity and they percolate at T^. 
The appearance of these superconducting droplets above is in agreement and it is the 
only possible explanation of various measurements made in the normal phase of different 
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46, 



47 



48|. In this 



materials like the Nernst effect 0, Q| and the precursor diamagnetism 
scenario, superconducting phase coherence is achieved only at Tc which is the temperature 
that ^ 60% (~ the percolating limit) of the sample volume is in the superconducting phase 
as has been proposed by several different works^, 0, |^. In the following sections we 
discuss the phase separation mechanism, we present the results of some simulation and the 
implications to HTSC properties in detail. 

As mentioned above, the process of phase separation in HTSC is well documented but, 
concerning the mechanism of phase separation there are not many conclusive studies. One 
possibility for this mechanism arises from the measurements by nuclear magnetic resonance 
(NMR)j2^, which has determined the high mobility of the oxygen interstitial in La2Cu04-s 
compounds. Therefore, it is possible that the dopant atoms cluster themselves to minimize 
the local energy an this would be a possible explanation for the whole process. This is just a 
general idea based on the NMR results j26| but the mechanism of clustering is an interesting 
subject that merits more attention in the future. 

To avoid confusion in the notation, we will adopt Tps{p) for the large pseudogap tem- 
perature of a compound with average hole doping p and T*{p) for the lower pseudogap 
temperature. When we refer to a given sample and not to a family of compounds, to sim- 
plify the notation, we may just use Tpg and T*. 
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II. THE CH APPROACH TO PHASE SEPARATION 



The CH theory was developed to the binary allows and one may question its application to 
a strongly correlated system as HTSC. However the clustering process in hole doped HTSC 
is very subtle. As we can draw from the stripe phases, the antiferromagnetic insulating 
phase has nearly zero holes per copper atom and the charged phase has less than 0.25 holes 
per copper atom, and in some cases 0.125. Thus, double hole occupancy does not occur in 
either phases, which is in agreement with a large on-site coulomb repulsion used in almost 
all Hamiltonian models for HTSC as in Eq.(@]) below. Therefore, we believe that the use of 
the CH theory to hole doped HTSC is justified. 

As an initial condition, let's suppose that a typical HTSC has, above Tpg, a Gaussian 
distribution of local densities around an average hole density p as can be direct inferred 



from the STM experiments 
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have measured a spread of Ap ^ 0.08 



33|. Pan et al 

holes/Cu for an optimally doped compound which will be adopted as an initial condition 
in our calculations. This Gaussian distribution around the average hole density p is the 
starting point at temperatures above and near the phase separation temperature Tpg and 
each local hole density p{x) inside the sample oscillates around the compound average p. In 
this way, we can define the order parameter u{x) = p{x) —p and u{x) = above and at Tps, 
as expected. Then the typical GL functional for the free energy density in terms of such 
order parameter is 

f = ^e'\Vu\' + V{u) (1) 

where the potential V{u) = A^{T)u^/2 + B^u^/A + A'^{T) = a{T - Tp,) and B is a 
constant. Notice that near and below Tp^ and/or for small values of e, the gradient term 
can be neglected and we get the two minima of / at the equilibrium values u{x) = ±A/B = 
±y^{a{Tps - T))/B. This can be easily seen if we write V{u) = B'^{u'^ - A^jB'^f. On 
Fig.(^ we show the important characteristics of such potential: As the temperatures go 
down away from Tp^, the two equilibrium order parameter (or densities) go further apart from 
one another and the energy barrier between the two equilibrium phases Eg also increases. 
Eg = A^(T)/B which is proportional to {Tpg — T)" 



Bray 4j| pointed out that one can explore the fact that the type of order parameter used 
above, as the two types of atoms of a given alloy, is conserved and the CH equation can be 
written in the form of a continuity equation, dtu = —V.J, with the current J = MV {6 f / 6u) , 
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FIG. 1: The typical potential used in the density of free energy which gives rise to phase separation 
as function of the order parameter u. Notice that the two minima at u± yield the two equilibrium 
densities p± = u±+p and the energy barrier between them Eg depend on the temperature difference 
T — T 

where M is the mobility or the transport coefficient. It is probably the same for each family 
of HTSC compounds because of the universal character of their phase diagram. Therefore 
we may write the CH equation as following, 

^ = -MV\e''V\ + A\T)u - B\'). (2) 

This equation is solved with the so-called fiux-conserving boundary conditions, 
Vu.n\^<zgn = (V^u).n\rg^dn = where n is the outward normal vector on the boundary 
of the domain Q which we represent by dQ, it is possible to show the time conservation 
of the total mass Mt and that the total free energy can only decrease (dissipate) or being 



stable |43|, l44l |. Therefore a time stepping finite difference scheme is defined to be gradient 
stable only if the free energy is non-increasing and gradient stability is regarded as the best 
stability criterion for finite difference numerical solutions of such non-linear partial differen- 
tial equation as the CH equation 14411. 

nn 

As it has already been pointed out|42l, |4^, both the and the non-linear term make 
the CH equation very stiff and it is difficult to solve it numerically. The non-linear term 
in principle, forbids the use of common Fast Fourier Transform (FFT) methods and brings 
the additional problem that the usual stability analysis like von Neumann criteria cannot be 



used. These difficulties make most of the finite difference schemes to use time steps of many 
order of magnitude smaller than Ax and consequently, it is numerical expensive to reach the 
time scales where the interesting dynamics occur. To solve these difficulties Eyre proposed 
a semi-implicit method in time that is unconditional gradient stable when the V{u) can be 
divided in two parts: V{u) = Vc{u) + Ve{u) where Vc is called contractive and Ve is called 
expansive l^]. Thus, we adopt here his method taking Ve as the quadratic term and Vc as 
the forth order one. Then we finally obtain the proposed finite difference scheme for the CH 
equation which is linearized in time (we have absorbed M into the time step), namely 

= U-,-/\tA\T)V'U^^,. (3) 

We have studied the stability conditions of this equation in one, two and three 
dimensions j^]. In the next section we present the results for two and three dimensions 
applied to the problem of phase separation in a HTSC plane of CuO. Although we calculate 
the local order parameter u{x) of a sample with average hole density p, we are interested 
and will preferably refer to the local hole density p{x) = u{x) + p. 



III. THE RESULTS OF THE SIMULATIONS 



As mentioned in the introduction, there is a consensus from several different 
experiments that the pseudogap temperature T*{p) initiate at average hole doping 
p ~ 0.05 at T ^ 800K and falls to zero temperature at a critical doping p^ ~ 0.2. This is 
best illustrated by Fig. (11) from the review work of Tallon et ai0| with many different data, 
which we reproduce here for convenience. 

Initially, that is above Tps{p), the system has a homogeneous distribution of charge with 
very small variations around p, which is described by a very narrow Gaussian-type distribu- 
tion. When the temperature goes down through Tps the sample with average hole density 
p starts to phase separate and the original Gaussian distribution of holes changes continu- 
ously into a bimodal type distribution. For underdoped samples with large Tpg, the mobility 
M is high which favors a rapid phase separation into two main hole densities p- and p+, 
while the compounds near the critical doping pc may not undergo a complete phase sepa- 
ration. Near the Tps, the difference between p_ and p+ is very small and increases as the 
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FIG. 2: Fig. 11 from Tallon et al;2'] showing the p-dependence of the pseudogap energy Eg or T* 
determined from susceptibihty, heat capacity, ARPES, ^^Y_NMR and resistivity as displayed in 
the legends. 



temperature goes away from Tpg. However if the system is quenched very rapidly the phase 
separation may even not occur, because it depends on the mobility which is essentially the 
phase separation time scalej4ll. [4^. For p > Pc there is no phase separation and the charge 
distribution remains a Gaussian like. For p < pc, the transformation from a homogeneous 
phase to one with different densities and with sites at different environment is seemed by 
many different measurements: by local measurements like, for instance, the Y3IMR, by 
transport measurements like the resistivity since the charges must overcome the potential 
barrier Eg between the two equilibrium regions (see Fig.((T))) and by susceptibility due to the 
appearance of antiferromagnetic regions with low hole density specially at the low average 
doping compounds. Notice that the coefficient A{T) = \/{a{Tps — T)) changes smoothly as 
the temperature goes down away from Tpg and therefore the charge distribution in a given 
compound depends strongly on the temperature T, on the details of sample synthesis and 
annealing procedures and, due to the mobility, on how the system is quenched through Tpg. 
This is probably the explanation to the different results reported in the literature on many 
HTSC compounds. ^ 

Assuming that the curve proposed by Tallon and Loram[2j reproduced here in the Fig.(j21) 
is the Tps line, the regions below are characterized by their temperature distance from this 
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temperature. The regions in the bottom hke 5 to 7, as illustrated in Fig.Q, are region with 
very strong phase separation while regions near Tps like 1 to 3 the phase separation is weak. 
This is because u± = ±iA/B) = ±^{a{Tps - T)/B) and these regions are characterized 
by their values of {Tps — T). Thus, in region one, the difference between p_ and p+ is very 
small and increases as the temperature goes below the Tps line. Accordingly, the energy 
gap Eg = Eg{T) is a varying function of T and goes to zero near Tpg. At zero temperature, 
compounds with p < 0.1 may be strongly separated in a insulator phase (p_ ~ 0) and in 
a metallic phase with > 0.2. Compounds with 0.1 < p < 0.16 the phase separation is 
partial and for 0.16 < p < 0.2 the original Gaussian is distorted with an increase in the hole 
density at the low and high tail. 



1000 
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FIG. 3: Illustration of the phase separation regions. The thick line represents Tpg or Eg from 
Fig.© approximated by a straight line. The numbered regions are equidistant from Tps and are 
characterized by their single values of (Tps — T) which is proportional to the equilibrium densities 
P- and p+. 

We have performed calculations in all regions below the phase separation line increasing 
the value of the A coefficient simulating the the temperature difference {Tps — T). Different 
initial conditions were tested to check convergence after thousands of time steps. One of the 
trial starting initial condition was, for instance, u{t = 0) = e x sin{x)sin{y). 

In Fig.(^ we show the results of the simulations on a 100 x 100 square grid. In these 
simulations we used A/ B = 0.125 and e = 0.05 which represents a phase separation in region 
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4 of Fig.Q because it is a region where phase separation is neither minimal as in region 
1 nor maximal as region 7. The simulation describes the time evolution of a homogeneous 
initial condition given above and represented by a very sharp Gaussian around the average 
p value shown in Fig.(j31|^a)). Fig.(j^ shows very clearly the phase separation process. 
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FIG. 4: The process of phase separation with the time. On the top panel we plot the order 
parameter map at the times t = 2000(a) which displays an (enhanced) reminiscent pattern from 
the initial conditions, At t = 5000(b) the phase separation process has already started and on the 
c) panel at the t = 10000 and t = 25000(d). 

The phase separation time evolution is also well illustrated by displaying the histogram 
of how the order parameter evolves in time. In the Fig.(0) we show the time evolution of a 
typical simulation with the same parameters of Fig.(0): t = 1 represents the initial condition 
with the hole density p{x) centered around an average value p = 0.125, t = 5000 represents 
5000 time steps in our simulations and so on. It is very interesting the shape of these 
histogram and their evolution from a initial centered Gaussian to a bimodal distribution. It 
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is very important to emphasize that the distribution after a certain time steps is independent 
of the initial condition. In practice, if the mobihty would be very large and if e is very small, 
the system would evolve to two delta functions at p±. 




FIG. 5: The evolution of the local densities of order parameter u{x) with the time in our simulations. 
We can see the tendency toward sharps bimodal distributions at the density equilibrium values 
{u- and uj^). 

To study the effect of the gradient term in the GL free energy of Eq.(0) we have also 
performed simulations with different values of e. We have tested e = 0.01, 0.03 and 0.05. 
The results are shown in Fig.® and we can see that indeed the order parameter distribution 
approaches a delta function as e decreases. 

Phase separation always occur when we start with a small variation around an average 
value but the final pattern is strongly dependent on the size of the system. In order to 
study such effect we have also done, together with the 100 x 100 lattice, calculation with the 
200 X 200 and 500 x 500 square grid. At the Fig.((7j), we show the results of mapping the order 
parameter in a surface with the same values of parameter used above. It is very interesting, 
in the context of HTSC, to observe that smaller lattices display granular pattern and there 
is a clear increase in the formation of stripes pattern as the size of the lattice is enlarged. 
It is a matter of fact that the largest HTSC single crystals are those of the Lai-xSrxCu02 
family which are more suitable for neutron diffraction studies and it is exactly in this family 
which the stripe phases were measuredjs^, 31 1. As conjectured by A.Moreo et alS^I, it is 
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FIG. 6: The evolution of the local densities of order parameter probability with the gradient 
constant £. We can see the tendency toward sharp bimodal formation at the values p± as e 
decreases. 

likely that the same conclusion may be applied to the manganites. 

Notice how the stripe structure develops in the plane interior and as they end at the 
borders they display a granular type pattern similar to those found in STM0, [s^ Issj ]. 
In order the check this we have also performed simulations in 3D. The results does not 
differ appreciably from the 2D case. At the Fig.® we show cuts in a three dimension 
100 X 100 X 100 lattice at the middle plane {z = 50) and near the top surface {z = 100). 
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IV. THE LOCAL GAP 



We have shown that that below the Tps a phase separation develops creating a variable 
density of holes at very small mesoscopic scale. Therefore, it is important to perform a local 
superconducting gap calculation, taking into account this charge inhomogeneity, in order to 
understand its effect on the the superconductivity phase and specifically, how such phase is 
built in this inhomogeneous environment. The appropriate way to do this calculation, in a 
system without spatial invariance, is through the Bogoliubov-deGennes (BdG) mean-field 
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FIG. 7: The mapping of the order parameter in the process of phase separation in lattices with 
different sizes. We display here the order parameter in the 200 x 200 (a) and for 500x500 (b) 
lattice. The parameter are the same and therefore it is to be compared with the results displayed 
in Fig.4d above for the 100 x 100 lattice. 



theory 
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57| . We start with the extended Hubbard Hamiltonian 



Ui^Uj^, (4) 



where cl^{cia) is the usual fermionic creation (annihilation) operators at site Xj, spin cr{ti}, 
and rijo- = cj^Cjo-- tij is the hopping between site i and j, U is the on-site and V is the nearest 
neighbor interaction, /i is the chemical potential and V^"^^ is a random potential which 
controls the strength of the disorder and introduces the inhomogeneous Hartree shi 
Using a mean-field decomposition approach, one can define the pairing amplitudes 
A5(xj) = V{ciiCi^s^) and A;7(xj) = [/(ci|Cj|), which yields an effective Hamiltonian 




i6a icr 



(5) 



In this expression 5 represents the nearest neighbor vectors and fii = fj. — j{ni) is the 
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FIG. 8: The 3D mapping process of phase separation at t=50000 steps .On the a) panel we plot 
the order parameter at the center z = 50 of a 100 x 100 x 100 lattice. At the b) panel we show the 
order parameter near the top surface {z = 100) 

Hartrcc shift with the local eletronic density (rii) = '^^{nia)- The hole density is p(xj) — 
1 — {rii). The H^ff is diagonalized by the BdG transformation 



En[7nT«n(Xi) -7li<(Xi)], 



(6) 



where 7^0. and 7^^. are quasiparticle operators associated with the excitation energies {E^ > 
0). Uni'X-i) and i'n(xj) are normahzed amphtudes for each Xj. Therefore the BdG equations 
are 



K A 



\ 



E„ 



(7) 



A* -K* 



\^^n(Xi) j 
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with 

KUn{:>^i) = - ^ti,i+5M„,(Xj + 6) + {V-"^'^ - jli)Un{y^i) 

5 

Au„(xi) = ^ A5(xi)u„(xi + (5) + A[7(xj)u„(xi), (8) 

and similar equations for f„(xi). These equations give the amphtudes ('U„(xi), t'„(xi)), and 
the eigenergies En- The BdG equations are solved self-consistently together with the pairing 



amplitude j5J, |55 1 



At/(xi) = -t/^M„(xiX(xi)tanh^^, (9) 



Asi^i) = -— J^[M„(x,X(xi + S)+ <(xi)M„(x, + 6)] tanh (10) 
and the hole density is given by 

P(X,) = 1-2 J2[Ki^^)\'fn + \Vni^i)\'{l ' fn)], (H) 

n 

where /„ is the Fermi function. Depending on the values of the potentials V and U, it is 
possible to have pairing amplitude with either s or wave symmetry 0, 0, 0|. 

It has been shown that a superconducting gap with d wave symmetry calculated in a 
square lattice, can be written as 

Ad(xi) = ^[Ax(xi) + A_5c(xi)-A^(xi)-A_^(xi)]. (12) 

Therefore we have used the above BdG theory to calculate the local superconducting zero 
temperature doping depended s and d wave gap. In Fig.Q, we show a typical set of results 
for d wave as function of doping, with y*"^? = and for a cluster of 14x 14 sites. We have used 
parameters which are appropriated to the HTSC, as we discussed in some of our previous 



works l5^ 



59[: a hopping value of t = 0.35eV, next neighbor hopping t2 = 0.55t, an on-site 
repulsion U = 1.3t and a next neighbor attraction V = — 1.6eV. Changing these parameters 
the gap curve also changes but its qualitative form is not affected. This calculation is to 
be used concomitantly with the phase separation results from previous sections, since below 
Tps, the system has regions or islands of different doping levels. The consequences of the 
BdG calculations, like those presented in Fig.Q, will be discussed in the next section in 
order to support the interpretation of many physical properties associated with the HTSC. 
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V. DISCUSSION 



As discussed in the introduction, it is very likely that phase separation is a fundamental 
process in the HTSC physics and therefore it must manifest itself through many experimental 
results. In order to explore this fact, we have developed a formalism based on the CH 
differential equation which allows one to quantitatively study the HTSC phase separation 
process. We take the upper pseudogap as the onset of phase separation because it starts 
in the underdoped region usually at very high temperatures (~ 800K) where we expect 
neither Cooper pair formation nor fluctuation of these pairs and also because there are many 
arguments against its identification with the superconducting gapQ]. In fact, the difficult 
to associate the experimental data at such high temperatures with superconductivity led 
some authors to call it simply a crossover temperature l|, Thus, assuming that the upper 
pseudogap temperature line as that shown in Fig.© is the onset of phase separation, we 
have been able to provide a simple interpretation to the occurrence of a gap (Eg) at such 
high temperatures, to follow the hole density time evolution and how a small fluctuating 
(almost homogeneous) phase separates into two main local densities (p_ and Now we 
want to discuss some more speciflc implications to the physics of HTSC if, in connection 
with the above, we take the lower pseudogap as the local onset of superconductivity. 

The lower pseudogap has been attributed to the local mean fleld (MF) superconduc ting 
temperature or to the onset of pair formation or superconducting fluctuation |6l. Il2l Il9l l5Cll| . 
Starting in the underdoped region at temperatures usually near the room temperature, 
it has been identifled with the onset of local superconductivity or with the appearance 



of small superconducting regions 5l|. This interpretation is supported by many different 
experiments, the most directs being the Nernst effect 0, Q| and muon spin rotation[19,j. 
Following the theoretical predictions |S] and the Nernst effect resultsj^, [3], we assume that 
the lower pseudogap vanishes at the strong overdoped region. Thus, in order to match the 
lower pseudogap, we have calculated the zero temperature superconducting gap A{p{x)) 
with a d-wave symmetry as described in the previous section and displayed in Fig.Q. It is 
important to use small clusters like 8 x 8, 12 x 12 and 14 x 14 in order to assure that we are 
indeed calculating the local properties but which are larger than the coherence length. 

It is interesting that the results of the local BdG zero temperature gap function A(x) 
have *e sa,ne ,ua..a«ve fonn of >owe. pseudogapQ B Q, and . y,eM. .a,e values 
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FIG. 9: Results of the calculation for the zero temperature local superconducting gap using a d- 
wave BdG superconducting theory. Since these calculations used a small clusters, we can attribute 
the average density to a local density p{x). 

at low doping with its maximum near p{x) ~ 0.05 and decreases continuously down to 
zero at the overdoped region. If the A{p) gap measured in a compound with average hole 
density p is assumed to be the corresponding average value of all A{p{x)), we arrive that the 
A{p{x)) xp[x) is very similar to the A{p) xp curve. Indeed the heat capacity measurements 
(see Fig. (8) of Tallon and Loram^) and the ARPES (see Fig. (4) of Harris et al[l^) yield 
A{p) X p curves with the same qualitative form of the A{p{x)) x p[x) shown in Fig.Q. 

Thus, the lower pseudogap temperature, which we denote T*{p), is the onset of supercon- 
ductivity and the superconducting regions grow as the temperature is decreased below the 
T*{p), but long range order is only possible at the percolation limit among these regions, 
when phase coherence is establish at Tc{p). This scenario, with these two (phase-separation 
and local superconducting) pseudogaps, is appropriate to interpret many non-convention 
HTSC features and their main phase diagram, that is, the curves Tpsip) (the upper pseudo- 
gap temperature), T*{p) (the (superconducting) lower pseudogap temperature) and Tc{p), 
as we discuss below: 

We start with the discussion of the many tunneling experiments results Q, Q, 
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one of the most well known fact about these experiments is that they do not yield 
and special signal at Tc{p) and form a kind of "dip" that persists above Tc{p) and dies off at 
T*{p). Our main point is that these experiments are made over a finite region and always 
measure the average of all A{p{x)) in this region. As the temperature is continuously raised 
from near zero, the regions with weaker A{p{x)) (and lower Tc{p{x))) become initially normal 
and, increasing more the temperature, many regions gradually turn from superconducting 
to normal state but, all the local superconducting regions are extinguished only at T*{p), 
not at Tcip). From the BdG calculations displayed in Fig.Q, we see that the regions with 
p near yield gaps near the minimum value A(p+) and we call it the lower or weaker 
branch. All the A(p(x)) in this branch which has their local densities p < p{x) < p+, 
vanishes before the temperature reaches Tc{p) while those in the strong branch A(p_) with 
P- < p{x) < p decreases also continuously as temperature is raised but they are more 
robust and totally vanishes only at T*. These features are probed by tunelling experiments 
which, due to the very small mesoscopic p{x) regions, usually measure the average of all 
these gaps. At low temperature, the average of many different gaps are measured and as 
the temperature is raised they all decrease and, firstly those in the weaker A{p+) branch 
and the ones in the second A(p_) afterwards, vanish at different temperatures from zero, 
passing by Tc{p) up to T*{p). Since all different A(p(x)) vary continuously, there is not any 
special or different signal at Tc{p). The measured dl/dV "dip" signal which came mostly 
from the robust gaps in the A(j9_) branch remains at temperatures well above Tc{p) in the 
underdoped region, decreases for compounds with increasing doping level p since p+ and p_ 
approaches one another, and in the overdoped region remains just for a few degrees above 
Tc{p). On the other hand, at the far overdoped region and above the critical doping p^ ~ 0.2, 
there is no phase separation and the distribution of p{x) is just a Gaussian-like distribution 
around the average hole doping p, consequently, the " dip" structure remains only for a few 
degrees, proportional to the distribution width. This is well illustrated by Fig. (3) of Suzuki 

n n 

et alp^ or in the figures of Renner et al[lJ|. As mentioned, these experiments see usually 
the average of many gaps in a given region, but, more recent refined experiments of Krasnov 
et alQl and Yurgens et al [23] were able to distinguish between the gaps in the two average 
A(p+) and A(p_) branches. Fig.(l) of Krasnov et al 22| shows clearly the weaker A(p+) (or 
superconducting in their interpretation) peak fading away as the temperature approaches 
Tc{p) while the larger average A(p_) (or their pseudogap) "dip" is almost unchanged. They 
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have also shown how apphed magnetic fields up to 14T destroy the weaker A(p+) leaving 
again the stronger A(p_) branch untouched. Thus, the pseudogap signal remains after the 
pair coherence is lost because the isolated or local superconducting regions left above T^. 
are those with very large A{p{x)) or Tc{p{x)) < T*{p) and the temperature and fields to 
destroy the superconductivity in these regions are much larger than Tc{p) and the 14T used 



in the experiment 



22| | . A similar finding was provided by a STM experiment which measured 



a remaining pseudogap signal inside cores of Bi-2212 quantized vortices, where long range 
superconductivity is clearly destroyedjc^. Furthermore, we predict that the average A(p_) 
maximum peak decreases slowly as the temperature tends to Tps because p+ and p_ coalesce 
to p. This temperature decreasing was recently measured and can also be seen in Figs. (2 
and 3) of Yurgens et al {2^ . 

These results and our interpretation also agrees with the high magnetic field experiments 
which have measured simultaneously the closing of the pseudogap field (Hpg) and T*{p) by 
interlayer tunneling and resistivity Their reported results are for compounds in the 
overdoped regime with p > Pc ~ 0.2, that is, for compounds with doping level above the 
phase separation critical doping and T*{p) is just the maximum local Tc{p{x)) or lower 
pseudogap. At these doping levels there is no phase separation and what Shibauchi et al 61 1 
measured as Hpg is the field that closes the maximum local critical temperature which is 
T*{p) because it is the largest of all locals Tc{p{x)). As they apply a magnetic field at low 
temperature, it destroys first the superconducting clusters with low local Tc{p{x)) and, as 
the field increases, regions with larger values of the Tc{p{x)) are destroyed. Increasing even 
more the external field, eventually it destroys the long range order or percolation among the 
superconducting regions at the superconducting close field Hgc, leaving still some isolated 
regions which have larger Tc{p{x)) than the phase coherence temperature Tc{p). Increasing 
more the field, one reaches the closing field Hpg = 60T which destroys all the superconducting 
regions at T*{p = 0.2). The closing field Hpg must be very high for compounds with lower 
doping, probably much higher than the 60T used b y S hibauchi et al|21| for a p = 0.2 
compound and that is the reason why Krasnov et alj22] did not see any change in their 
optimally doped pseudogap dip at 14T. On the other hand, a d-wave BCS with a Zeeman 
coupling yields good agreement with the data, supporting the origin of the lower T* as 
the maximum local superconducting T(.{p{x)) 2u\. The fact that the local superconducting 
regions with large Tc{p{x)) and low local doping (around p_) are very robust to an external 
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magnetic field is also consistent with the Knight shift measurements which have seen the 
reductions of 1/TiT and K above from the values expected from the normal state at high 
temperatures in the overdoped region without any field effect up to 23. 2T in the underdoped 
region 6^. Notice that, since the experiment of Shibauchi et aljo^ is performed with p > 
0.2 samples, that is above the phase separation threshold Pc, therefore there is only one 
(Gaussian) dispersion of local superconducting gaps A(x) and there is no gap Eg associated 
with any phase separation. 



More recently, Hoffman et al|63l |6J and McElroy et alj65|, |66l developed a refined STM 
analyses which let them to study the doping dependence and the electronic structure of some 
compounds of the Bi-2212 family jc^. They find a distribution of low temperature local 
superconducting gap values A{x) whose average value A(p) and its width at half maximum 
increases for compounds with average hole doping varying between p ^ 0.19 and p ~ 0.11. 
The measured local values of A(x) varies from 20meV to 70meV at regions with linear sides 
of approximate 55nm in length. The low energy gaps exhibit periodic modulations consistent 
with charge modulations like of a granular charge phase separation. Their results, specially 
those shown in their Fig.(3A-3E), display a distribution of mesoscopic scale regions local 
gaps of two types: 

a) One type derived from a dl/dV curve with sharp edges with values < 65meV which 
they called coherence peaks. They interpreted this type of peaks as due to superconducting 
pairing on the whole Fermi Surface arguing that this kind of spectrum is consistent with a 
d-wave superconducting gap 

b) Another type of dl/dV spectra display an ill defined edges of a V-shape gap with larger 
values than +65meV what they called zero temperature pseudogap spectrum. Furthermore, 
they find that the a-type spectra are dominant for overdoped samples {p = 0.19 and 0.18) 
in which there is practically 0% probability of occurring spectra of b-type. The b-type 
spectra start to have a nonzero probability for compounds with p = 0.14 or below, and for 
underdoped compounds like p = 0.11 they find more than 55% of b-type spectra. It is not 
difficult to explain these observations in terms of the CH phase separation scenario: The p = 
0.19 and 0.18 compounds are near the phase separation threshold and their p{x) distribution 
is essentially a Gaussian type. Eg is small and they measure a Gaussian distribution of local 
superconducting gap values A{x) or a-type spectra. On the other limit, for the p = 0.11 
compound, according to the phase separation histogram of Fig.©, almost half of the system 
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has very low doping (0 < p{x) < 0.5), the j9_ branch and almost half of the system is 
in the other p+ branch (0.17 < p{x) < 0.22). The regions with p{x) in the p+ branch 
exhibit a superconducting gap distribution of a-type spectra, while regions in the p_ branch 
are mainly in the insulating region which produces b-type spectra. These a- and b-type 
spectra are mixed in the intermediated p = 0.15 and p = 0.13 compounds and it is in this 
near optimal doping region that both superconducting and the b-t ype gap have equally 



proba 



al 



65 



3ility. Consequently, the gap maps found by Hoffman et alj63|, |6J| and McElroy et 



66l | are a clear manifestation of the phase separation process in Bi2212. 



There are many other experiments which we could discuss in the light of the present phase 
separation theory but we believe that the above discussion is sufficient to demonstrate that 
phase separation process is central to understand many non-conventional HTSC properties. 



VI. CONCLUSIONS 



We have studied analytically the problem of phase separation in HTSC taking some cur- 
rent ideas on the possibility to identify the upper T* with the onset of phase separation 
and the lower pseudogap as the onset of d-wave superconductivity. Our approach allows us 
to make quantitative calculations of the phases separation process and to perform simula- 
tions which led to granular and stripe patterns depending on the parameter and size of the 
lattice, which are in agreement with current observations. Such calculations might be also 
be pertinent to the physics of manganites. It is also possible to get some insights on many 
general experimental results like: 

i- The charge distribution becomes more inhomogeneous in the underdoped region of the 
phase diagram where the stripes has been observed because Tp^ is larger in this region. 

ii- The spatial variation or width of the local hole concentration p{x) increases as the tem- 
perature decreases. 

iii- The fact that some materials exhibit granular while others exhibit stripe patterns may 
be related with the single crystal or palette size of ceramic or granular samples. Our sim- 
ulations indicate that larger lattices favor stripe patterns while smaller ones favor granular 
patterns. 

iv- The spinodal decomposition reveals the importance of the sample preparation process, 
that is, samples with the same doping level may have different degree of inhomogeneity 
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depending on the way they have been quenched through Tpg. This would explain different 
results on the same kind of compounds which has been very frequent in the HTSC 

v- The two different signal detected by refined tunelling experiments. 

vi- The density os state modulation measured by recent STM data. 

In summary, the CH phase separation approach to HTSC in connection with local charge 
density dependent BdG superconducting critical temperature calculation is used to explain 
the existence and nature of the two different pseudogaps and it provides novel interpretations 
on many non-conventional features and inhomogeneous patterns. Therefore, our main point 
is that we should regard the phase separation process as one of the key ingredients of the 
HTSC physics. 
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